Comparative study of defect energetics in Hf02 and SiC>2 
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We perform ab initio calculations, based on density functional theory, of substitutional and va- 
cancy defects in the monoclinic hafnium oxide (m-Hf02) and a-quartz (Si02). The neutral oxygen 
vacancies and substitutional Si and Hf defects in Hf02 and Si02, respectively, are investigated. Our 
calculations show that, for a large range of Hf chemical potential, Si substitutional defects are most 
likely to form in Hf02, leading to the formation of a silicate layer at the HfCb/Si interface. We also 
find that it is energetically more favorable to form oxygen vacancies in SiC>2 than in Hf02, which 
implies that oxygen deficient HfCh grown on top of Si02 will consume oxygen from the SiC>2. 
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The continuous device miniaturization in the mi- 
croelectronic industry will eventually lead, within the 
present technology, to the end of the use of amor- 
phous SiC>2 (a-Si02) as gate dielectric in metal-oxide- 
semiconductor field-effect transistors (MOSFETs). The 
existence of a thickness limit for the a-SiC>2 around 10-12 
A, has clearly been established experimentally pj. One 
way to circumvent this problem, still keeping Si as the 
basic device material, is to employ high-permitivity ma- 
terials as alternative gate dielectrics in place of the con- 
ventional a-SiC>2. Among them, hafnium oxide is emerg- 
ing as the material with greatest potential to substitute 
SiC>2, mainly due to its high dieletric constant and ther- 
modynamic stability, when it forms interface with Si. 

Even though hafnium oxide is thermodynamically sta- 
ble against an overall decomposition as Hf and SiC>2 when 
grown on Si, interfacial reactions can occur. Thereby the 
formation of a thin interfacial layer (oxides, silicates and 
silicides) between the HTO2 and the Si surface, has been 
recently observed I^Q. This interfacial layer occurs dur- 
ing almost any film growth processes or post-annealing, 
which is an intrinsic part of any growth cycle. Therefore, 
the thermodynamic stability of the hafnium oxide in con- 
tact with silicon is identified as a critical issue for the 
application of alternative gate dielectric in silicon-based 
devices 0, Moreover, the study of possible defects 
related to the migration of atoms across the interface is 
of fundamental importance. In particular, a significant 
source of defects in this system is the interface itself, 
which has been shown |(j to consist of Hf silicates with 
dielectric constant lower than that of HfC>2 0, |1| . 

In the present work we address the formation of neu- 
tral defects through first-principles calculations, based on 
the density functional theory (DFT). We analyze the for- 
mation of Si substitutional defects in HfC>2, as well as Hf 
substitutional defects in SiC>2, for different growth con- 
ditions. Finally, the energetics of an oxygen vacancy in 
SiC>2 is compared to a similar vacancy in HfC>2, in or- 
der to understand the growth of hafnium oxide under 
oxygen-poor conditions. 
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Many experimental works have addressed 

the chemical reactions that could occur in the HfC^/Si 
interface during the HfC>2 growth cycle. This is due in 
order to prevent and/or control the interfacial layer for- 
mation. Almost all works have reported the formation of 
an interfacial Hf silicate in oxygen-rich atmospheres. Fur- 
thermore, Wang and co-workers have shown that in 
opposite conditions, i.e., oxygen-deficient atmospheres, 
the Hf silicate interfacial formation does not happen dur- 
ing the HfC>2 growth cycle. However, there are still many 
open questions, such as which atomic species are migrat- 
ing when the interfacial silicate is formed. For example, 
it is important to know if Hf will be incorporated in a 
formed SiC>2 layer, or if Si from either the bulk Si or this 
SiC>2 layer will be incorporated in the Hf02- 

The DFT calculations were performed using ultra- 
soft Vanderbilt pseudopotentials 0], and the gener- 
alized gradient approximation(GGA) for the exchange- 
correlation potential as implemented in VASP code 
111 13 IH 111 In order to study the defects in the dif- 
ferent systems we have considered the monoclinic HfC>2 
and the a-quartz (SiO^) crystalline phases, using a 96 
atoms and 72 atoms supercells, respectively. For these 
cells, we have used a plane wave cutoff energy of 400 
eV and a 2 x 2 x 2 Monkhorst-Pack k-mesh. These crystal 
structures have been previously used to describe these 
systems 0, IH C01 ■ In a U calculations the atoms were 
allowed to relax until all components of the atomic forces 
were smaller than 0.025 eV/A. 

The silicon substitutional defect (SiHf) was created in 
the m-Hf02 supercell by substitution of one hafnium 
atom by one silicon atom in the equilibrium perfect crys- 
tal. On the other hand, the hafnium substitutional de- 
fect (Hfgi) was created in the a-quartz by substitution of 
one silicon atom by one hafnium atom in the equilibrium 
perfect crystal. Our results show that the presence of 
either Sinf or Hfsi in Hf02 or Si02, respectively, do not 
introduce any additional active levels in the band gap. 
For the perfect Hf02 at the equilibrium, we obtain an 
indirect band gap of 3.9 eV along T-B. For this same su- 
percell, the density of states (DOS) shows that the 0(2s) 
and 0(2p) bands are centered at around -14 and eV, 
with bandwidths of approximately 4 eV and 8 eV, re- 
spectively, below valence-band edge, whereas the Hi(5d) 



band is centered at around 8 eV forms in the conduction 
band, with bandwidth of 6 eV. These results are in good 
agreement with previous DFT calculations of this ma- 
terial [23. Moreover, for HfC>2 with a Si substitutional 
defect, we observe a resonant state around 6.9 eV be- 
low the top of the valence band, which is related to the 
0(^p)-like and Si bonding state |23| . 

The formation energy for a Sinf is calculated as 

Ef Hf = \E t (Si m ) +HHf]~ [E t {Hf0 2 ) + fi 8i ], (1) 

whereas for Hfsi the similar expression is 

Ef fs > = [E t (Hf Si ) + ^} ~ [E t (Si0 2 ) +HHf]- (2) 

In the above expressions, E t (D,g) are the total energies of 
the fully relaxed supercells (either m-HfC>2 or Si02) with 
the substitutional defect Dg, and Fi t (X0 2 ) are the total 
energies of the similar supercells for the perfect crystals 
of XO2 (X = Hf or Si). The values of the chemical po- 
tentials, fiHf and ^sij depend on the growth conditions. 
We have considered two limits for i^nf- (i) The bulk 
metal as a reference, which would correspond to a Hf- 
rich growth condition and the formation of Hf clusters in 
the bulk or at the surface of the oxide; (ii) Under oxy- 
gen rich conditions, or considering that there is plenty of 
oxygen atoms in the HfC>2 [24|], such that the removed 
hafnium remains always in equilibrium with the oxygen 
gaseous, the /x_y/(gas) can be obtained as /xjy/(gas) = 
^Hf0 2 ~ A*cv The HfC>2 chemical potential (a*h/0 2 ) was 
obtained as the energy per unit formula for the mono- 
clinic bulk hafnium oxide, and fio 2 is the energy of an 
isolated oxygen molecule, which was obtained through a 
DFT total energy calculation for an O2 inside a cubic 
supercell of 15 A side. For the silicon chemical poten- 
tial, two similar limits were considered: (i) /isi as the 
crystalline bulk Si chemical potential, and (ii) fjtsi as 
the chemical potential of Si in SiC>2 under oxygen rich 
conditions, i.e., fisi = MSi0 2 — Mo 2 ; with /io 2 as above 
and /J,si0 2 obtained as the energy per unit formula of 
a-quartz. 

The results for the formation energies are presented in 
Fig. ^ As can be seen from Fig. ^a), under oxygen rich 
conditions, the formation of a Si substitutional defect in 
HfC>2 is very likely to occur, specially if the Si source is 
the bulk crystal. Under this O-rich situation, the forma- 
tion of this defect would only become unfavorable if a 
thick enough layer of SiC>2 would exists between the Si 
substrate and the HfC>2. In this Case, either the source 
of Si atoms would be the silica, or the Si would somehow 
not diffuse through this thick layer. If the SiC>2 is thin 
enough, such that the Si diffusion through it is signifi- 
cant, then the formation of a Hf silicate (at least close 
to the interface) seems unavoidable [25|. On the other 
hand, under Hf rich conditions, it is less likely that Si 
substitutional defects in Hf02 will be formed. 

The formation of Hf substitutional defects in SiC>2 has 
a behavior that is opposite to Sinf , as can be seen in Fig. 
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FIG. 1: Formation energies for the two substitutional defects 
considered in this work: (a) Si in place of a Hf in Hf02, and 
(b) Hf in place of a Si in Si02- The formation energies are 
plotted as a function of the Hf chemical potential, and for 
two values of the Si chemical potential, the bulk Si chemical 
potential (solid curves) and the chemical potential for Si in 
the Si02 under an oxygen rich environment (dashed curves). 

1(b), i.e., it is not likely to form under oxygen rich con- 
ditions, and becomes more probable under Hf-rich con- 
ditions. In a situation where the Si chemical potential is 
given by its bulk value (either thin or no Si02 layer), the 
formation of a Hfgi in Si02 is unlikely to happen, for a 
large range of Hf chemical potential. 

We also considered the formation of neutral oxygen 
vacancies (Vo), both in SiC-2 as well as in HfC>2 [2j], since 
they can be created in films and bulk samples due to the 
growth cycle. The neutral oxygen vacancy in the ?n-Hf02 
and a-quartz were generated by simple removal of an 
oxygen atom, followed by full relaxation of all remaining 
atoms. The formation energies for a Vo in XO2 (X=Hf 
or Si), Ef(Vo), were calculated as, 

Ef°*{V ) = [E?°*(V ) + no] [E t {X0 2 % (3) 

where E t (Vo) and E t (X0 2 ), are the total energies of 
supercells of X0 2 (X = Hf or Si) with, and without, an 
oxygen vacancy, respectively. The oxygen chemical po- 
tential, /io, was considered either as the total energy of 
an isolated oxygen atom, or as one half of the energy of 
an isolated oxygen molecule (in both cases, the O2 triplet 
ground state was used). The monoclinic phase has non- 
equivalent oxygen atoms, i.e, which some sites are three- 
fold coordinated, whereas others are fourfold coordinated 
by hafnium atoms. In this way, we have determined the 
formation energy for both vacancies types using Eq. © 
and we obtain a formation energy difference around 0.02 
eV. In the former case, we obtained of Ef f ° 2 (V ) = 
9.32 eV and Ef° 2 (V ) = 8.10 eV, whereas for the lat- 
ter choice of ^o, we obtained E^° 2 (Vo) = 6.38 eV 
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and E s f l ° 2 (V ) = 5.16 eV. This indicates that, although 
neutral oxygen vacancies are energetically unfavorable in 
both materials, they are more stable in silicon oxide than 
in hafnium oxide, by approximately 1.23 eV. This implies 
that, if an oxygen deficient Hf02 is grown on top of a SiC>2 
layer, oxygen atoms will migrate from Si02- In this way, 
oxygen vacancies are created in the silica layer toward the 
HfC>2, healing at the sometimes oxygen vacancies in the 
hafnia. Indications that this process does indeed occur, 
has been recently reported [l3j . 

In summary, our results show that, unless the hafnium 
chemical potential is always very close to its bulk 
value, i.e., oxygen-poor growth conditions, the formation 
of Si substitutional defects in the HfC>2 is almost unavoid- 
able. This will lead to the formation of a silicate-like layer 
close to the Interface. For a very large range of Hf chem- 
ical potential it is very likely that incoming silicon atoms 



from the Si substrate will form a (Hf02)z(Si02), / layer, 
as supported by experimental observations 25] . The only 
way to prevent the formation of a silicate layer would be 
through the use of an oxygen deficient HfC>2, as has been 
observed [13J. However, even if this silicate formation 
is prevented in the initial steps of dielectric growth, it 
seems almost impossible that it will not be formed given 
the necessary further thermal processing steps, unless a 
barrier for the Si diffusion being introduced [9j]. 
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